The Evolution of the Intergalactic Medium Transmission to 

Redshift Six 



Antoinette Songaila^ 

Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822 

Received ; accepted 

To be published in Astronomical Journal, May 2004 



^Visiting astronomer, W. M. Keck Observatory, jointly operated by the California Insti- 
tute of Technology and the University of California. 



- 2 - 
ABSTRACT 

We have measured the transmission of the Lyct forest produced by neutral 
hydrogen scattering in the intergalactic medium between redshifts 2 and 6.3 using 
high signal to noise, high resolution {R > 5000) observations of 50 quasars spread 
over the redshift range. We use a uniform set of 15 A intervals covering Lya, 
Ly/3, and Lyj absorption regions to tabulate the forest transmission as a function 
of redshift. The transmitted fractions show a relatively smooth evolution over 
the entire range of redshifts, which can be modelled with a smoothly decreasing 
ionization rate. Previous claims of an abrupt change at 2; ~ 6 appear in part 
to be a consequence of an incorrect conversion of Ly/? to Lya optical depths. 
The tabulated transmissions can be used to calculate the colors of objects with 
a specified input spectrum as a function of redshift. We calculate the colors of a 
flat fi, galaxy with a large intrinsic continuum break, as an important example. 

Subject headings: early universe — intergalactic medium — quasars: absorption 
lines — galaxies: formation 
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1. Introduction 

The epoch of hydrogen reionization is one of the landmarks of the high redshift uni- 
verse, providing constraints on the first hght sources to have formed and a crucial input 
to models of structure formation. Recent dramatic progress has been made in determining 
limits on this redshift, leading to an interesting observational situation. On the one hand, 
spectroscopy of the Lyman a. forest shows that the intergalactic medium (IGM) is highly 
ionized at 2; < 6, becoming progressively opaque to radiation blueward of the quasars' Lya 
emission, culminating in the spectrum of the z — 6.28 QSO SDSS 1030+0524 which exhibits 
a substantial dark region that is consistent with no transmitted flux, i.e. a Gunn-Peterson 
trough at ^ = 6.05 (Becker et al. 2002, Songaila & Cowie 2002; Fan et al. 2002). This high 
opacity has been used to infer that this redshift is at the tail end of the hydrogen reionization 
epoch (Becker et al. 2002; Djorgovski et al. 2002) and comparison with simulations suggest 
the IGM could have been fully neutral at z = 6.5 (Cen & McDonald 2002; Gnedin 2002; Fan 
et al. 2002) though, as we point out in this paper, some of the strongest claims in this regard 
are based on an incorrect conversion of Ly/? opacities to equivalent Lya opacities (Becker et 
al. 2002; White et al. 2003). On the other hand, the recent detection by WMAP of a large 
optical depth to electron scattering (Bennett et al. 2003) is consistent with the Universe 
having been reionized at ;2 ~ 15 and fully ionized thereafter. 

The exact significance of the Gunn-Peterson trough in SDSS 1030+0524 is still not clear. 
It could follow from the ongoing thickening of the Lya forest rather than representing the 
onset of reionization since it is extremely difficult to distinguish the two effects. The discovery 
of three new 2; > 6 SDSS quasars (Fan et al. 2003) has made the situation more complex 
since, as might have been expected, there is significant variance in opacity between different 
hues of sight at ^ ~ 6 (Fan et al. 2003; White et al. 2003). Although White et al. (2003) have 
argued that apparent increased transmission in the spectrum of SDSS 1148+5251, currently 
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the highest redshift quasar, could be a result of foreground emission, it is more likely to be 
real cosmic variance (Songaila 2004a). In assessing the significance of the opacity at 2; = 6 
we need to know how the ionization rate in the IGM is evolving at lower redshifts, since the 
strongest signature of reionization would be an abrupt change in ionization at some redshift. 
In this paper, we use high signal-to-noise spectra of a large sample of quasars, including 
all the brightest z > 4.5 quasars observable from the northern hemisphere, to investigate 
the evolution of the transmission in the IGM from z — 2 to z — 6.3. We find that the 
transmission is smoothly evolving with redshift throughout the range. 



2. Observations 

Spectra of 25 high redshift quasars, with 4.17 < Zem < 6.39, were obtained for this 
program using ESI (Sheinis et al. 2000) on the Keckll telescope in echellette mode. The 
resolution was ~ 5300 for the 0.75" slit width used and the wavelength coverage is complete 
from 4000 A to 10, 000 A. Sample spectra and a more extensive discussion can be found 
in Songaila & Cowie (2002) and a complete catalog will be published in Songaila (2004b). 
At the red wavelengths, extremely high precision sky subtraction is required since the sky 
lines can be more than two orders of magnitude brighter than the quasars. In order to deal 
with this issue, individual half-hour exposures were made, stepped along the slit, and the 
median was used to provide the primary sky subtraction. Total exposure times were around 
3-7 hours per quasar, and are summarized in Table la. The frames were then registered, 
filtered to remove cosmic rays and artifacts, and then added. At this point a second linear 
fit to the slit profile in the vicinity of the quasar was used to remove any small residual sky. 
The observations were made at the parallactic angle and flux calibrated using observations 
of white dwarf standards scattered through the night. These standards were also used to 
remove telluric absorption features in the spectra, including the various atmospheric bands. 
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The final extractions of the spectra were made using a weighting provided by the profile of 
the white dwarf standards. Measurements of residual fiux in saturated regions of the spectra 
suggest that the sky subtraction leaves about 1% residuals in the faintest quasars and smaller 
errors in the brightest ones. We take the larger value as our systematic error. Wavelength 
calibrations were obtained using third-order fits to CuAr and Xe lamp observations at the 
begining and end of each night, and the absolute wavelength scale was then obtained by 
registering these wavelength solutions to the night sky lines. All wavelengths and redshifts 
are given in the vacuum heliocentric frame. 

A low redshift comparison was provided by spectra of 25 quasars obtained for other 
programs. This sample is summarized in Table lb. These spectra were taken with the 
HIRES spectrograph (Vogt et al. 1994) on the Keck I telescope with slit width 1.1 arcsec, 
giving a resolution R = 36, 000. Details of the data reduction have been described previously 
(Songaila 1998). 

The simplest characterization of the effects of the scattering by the forest is its averaged 
effect in depressing the continuum. This is usually characterized by the Oke indices (Oke & 
Korycansky 1982) averaged over the total Lya, Ly/3, etc. regions of the spectrum, or by the 
generalization to the transmitted fraction. 



where f\ is the observed spectrum, /cont is the continuum level in the absence of scattering, 
and the brackets indicate averaging over a specified wavelength region. 

Two methods have been used to compute /cont- For low resolution spectra, or for high 
redshift spectra for which F is very small, we can calculate /cont only by extrapolating from 
emission line-free regions in the spectrum longward of the quasar's Lya emission. Ideally 
one fits both the normalization and the power law index but in the highest redshift quasars 
there is seldom enough line-free spectrum to justify the fitting of the power law index, and 




(1) 
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we assume instead a fixed power law index (here —1.25) and set the normahzation using the 
region of the spectrum at a rest wavelength of 1350 A. When F is small, as is the case at 
high redshifts, it is relatively insensitive to the choice of power law extrapolation. At lower 
redshifts, where F is close to 1, the measurement is sensitive to the extrapolation, and here 
it is preferable to make a local interpolation between the peaks in the spectrum to determine 
the continuum (Zuo & Lu 1993). This in turn depends on there being regions of low opacity 
where the spectrum is close to the continuum value. 

We have used the continuum interpolation procedure for all of the HIRES data and 
for the ESI data at redshift less than z — 4.5, and have used the extrapolated continuum 
procedure for the ESI data at redshifts above 4.5. We have compared the two methods in 
the redshift range 4.0 < z < 4.5 where there are still enough low opacity regions to allow 
us to use the continuum interpolation method but where there is substantial absorption 
through most of the spectrum, so that the transmission is only weakly dependent on the 
power law extrapolation. For the 85 measurements in this range in the ESI spectra, we 
find that, for a power law extrapolation with index 1.25, the ratio of the measured value 
relative to the interpolated continuum transmission is 0.84 with a dispersion of ±0.18. For 
an extrapolation with a power law index of 1.75, this becomes 0.92 ± 0.20. The measured 
transmission fractions lie in the range 0.2 to 0.6 so this is a relatively small correction and the 
slightly lower value of the continuum interpolation method is caused in part by its tendency 
to over-estimate the transmission slightly because the continuum fitting regions have some 
opacity themselves. This effect becomes larger with increasing redshift. However, because 
the local continuum fitting method is more robust for individual quasars, wc have chosen to 
use it for the ESI data below z — 4.5 and it is these values that are quoted in Table 2. 

Fluxes were measured in regions of the spectrum away from the emission lines, namely: 
for Lya, in 8 bins of width 15 A between 1080 A and 1185 A; for Ly/?, in two 15 A bins from 
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980 A to 1010 A; for Ly7, in 1 bin of width 15 A from 955 A to 970 A. The mean transmitted 
fractions for Lya and Ly/3 are plotted in Figs, la and lb (small filled squares) as a function 
of redshift for the z > 4.5 ESI sample. Large open squares indicate quasars with BALQSO 
signatures, but the transmission of these objects do not appear to differ significantly at these 
wavelengths and we do not subequently distinguish these. The transmissions as a function 
of redshift are summarized in Table 2. The transmitted fluxes blueward of Lya and Ly/? in 
the HIRES sample, measured in the same rest-frame redshift bins as for the ESI sample, are 
plotted as crosses in Figures la and lb and tabulated in Table 2. Diamonds denote the ESI 
measurements with A < z < 4.5, made with the local interpolation method. Values of the 
transmitted flux in the redshift range z < 5 are in general agreement with the observations 
of Schneider, Schmidt & Gunn (1991) and Kennefick, Dgorgovski & de Carvalho (1995) and 
with the earlier hterature quoted in these papers. See also Songaila & Cowie (2002) for a 
comparison. At redshifts z > 5, they are also in general agreement with the results of White 
et al (2003). The Lya optical depths from White et al. are shown in Fig. 2a for the quasars 
SDSS 1148+5251 (open diamonds) and SDSS 1030+0524 (open squares). 

For purposes of display, it is useful to define a quantity tmean = — log(transmitted 
fraction)but it is important to remember that this is not the average of the optical depth. In 
particular, tmean for Ly/5 is not simply related to tmean ioi Lya but is a complex function 
of the structure (Fan et al. 2002; Songaila & Cowie 2002). Incorrect interpretation of this 
quantity as a simple optical depth, as has been done in Becker et al. (2002) and White 
et al. (2003), results in large overestimates of the Lya optical depth constraints based on 
Ly/5 transmission. Using the mean of the Lya and Ljf] transmissions, averaged in six-point 
bins, we have computed the corresponding Lya and Ly/? tmean as a function of the median 
redshift of each bin. These are shown in Figures 2a and 2b, where the error bars show the 
range of optical depths corresponding to the maximum and minimum values of transmitted 
flux in each bin. 
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3. Discussion 

A uniform ionized IGM would have a Lya optical depth 

65 km s-i Mpc-i J \ 7 J ^ \ 7 J 

(2) 

where r_i2 is the local ionization rate produced by the metagalactic ionizing flux in units of 

10~^^ s~^, T4 is the gas temperature in units of 10^ K, and the second part of the equation 
defines g, which is the normalized ionization rate. For such a medium, is simply = 0.16Ta, 
where the numerical factor is the ratio of (/A) for the two lines, where / is the oscillator 
strength. (Note that Becker et al. (2002), Fan et al. (2002) and White et al. (2003) incorrectly 
used only the ratio of the oscillator strengths.) The dashed lines in Figs. 2a and 2b show this 
shape for a constant F and T4 (or constant g), where we have chosen the normalization to 
provide a rough match to the Lya lines. This greatly oversimplified model gives a surprisingly 
good representation of the Lya transmission over the full redshift range but underpredicts 
Ly/3, as can be seen in Figure 2b. 

At 2; ~ 4.5 the average-density IGM begins to have substantial optical depth and it is 
at this point that very few regions return close to the extrapolated power law continuum 
as even the underdense regions have some significant opacity. It is also at this point that 
the transmission starts to fall extremely rapidly as the underdense voids themselves become 
opaque and there are no longer any regions of the IGM that transmit significantly. In this 
limit, it is clear why the ionized Gunn-Peterson formula is a relatively good approximation, 
since the bulk of the volume of the IGM is opaque and also relatively slowly evolving in 
filling factor and density. Because the transmission at high redshifts is dominated by the 
lowest density "voids" corresponding to the most underdense regions that occupy a significant 
volume, and since these regions evolve in a very slow and straightforward way, it is possible 
to characterize very simply the evolution of the mean transmission. Songaila & Cowie (2002) 
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derived the formula 

F{z,g) = 4.5^-°-^« (^)''' (-4-4 (^)') 

for the mean Lya transmission as a function of redshift for the ACDM model of Cen & 

McDonald. This analytic result closely matches the exact calulations of of McDonald and 
collaborators (McDonald & Miralda-Escude 2001, Cen & McDonald, 2002). It can be gen- 
eralized to calculate the transmission at Ly/3: 

Fp{z,g)^F{g,zp)F{g/0.16,z) (4) 

which is simply the product of the Ly/? transmission and Lya transmission at the redshift 

1 + Zf3 = {\/3/\a){l + z). Since we do not expect any substantial dependence on cosmology for 
any model that matches the lower redshift Lya forest, these equations allow us to express the 
evolution of the transmission in terms of the redshift evolution of the ionization parameter, 
g. For Lya, we use Equation (1) to fit the transmitted flux, assuming is a power law 
function of (1 + z), and obtain g = 0.74 x ((1 + z)/6)^'^'^. Following the prescription given 
above, the Ly/? transmission is then determined and is given by equations (3) and (4). For 
this more realistic model, the ratio of the Ly/3 and Ly a transmissions is 0.4, which is a 
factor of 2.5 times larger than the value of 0.16 which is valid for the ionized Gunn-Peterson 
case. As a consequence, models that fit Lya, such as the simple power law ionization we 
have assumed here, also provide a reasonable fit to Ly/3, as can be seen in Figs. 2a and 2b. 

However, rather than fitting simple models for the ionization, it is better to simply 
invert the equation (McDonald et al. (2000)), and in Figure 3 we have shown the ionization 
parameter at each redshift which would be required to produce the observed mean transmis- 
sion. Using Equation (3), we have calculated the ionization rate corresponding to the mean 
Lya transmissions computed in 12-point bins, starting with the highest redshift. These are 
displayed as the small filled squares in Figure 3. The horizontal error bars show the range of 
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redshift in each bin and the vertical error bars show the ionization parameter corresponding 
to ±1 a in the mean transmission in each bin. The open squares are the measurements 
of McDonald et al. (2000) at lower redshift. The dashed line is the power law fit to the 
presently measured values of for 2; > 4 given above. The simple formula of Equation (2) is 
valid only at z > 4, so we show the interpretation of the data only above that point. As long 
as Ly a transmission is detectable, the Ly a measurement should give the most sensitive 
and reliable measurement of g. However where Lyct is not seen, the Ly/? can give a more 
sensitive constraint because of its weaker oscillator strength (Becker et al. 2002; Fan et al. 
2002; White et al. 2003). We note that there is some sign in Fig. 2b that the Ly/5 opacity 
may be high relative to Lya in the highest redshift bin. The value of g would be 0.03 at 
z — 6.05 if we used Ly/3, which is more consistent with Fan et al. (2002) though it still lies 
above their upper limit. However, we emphasise that the Lya values are more robust. This 
drop off in the ionization parameter must also be combined with the mean free path for 
ionizing photons to determine the actual production rate of ionizing photons. This in turn 
means that the fall off in the rate of ionizing photon production is shallower than the fall-off 
in g. We postpone a more detailed discussion of this to a subsequent paper. 



4. High Redshift Colors 

Precise measurements of the forest transmission are also invaluable for determining the 
colors of high redshift galaxies and quasars. Currently the effects are usually included using 
the Madau (1995) models which are based on early data in the range 2 < 2; < 5. The present 
data can be used instead to compute these colors on an empirical basis and also to compute 
how much change can be produced by cosmic variance in the forest scattering. For broad 
band filters this latter effect is very small at all redshifts but it can be significant if narrow 
band filters are used. 
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As a useful example, we have computed the expected magnitude change as a function 
of redshift for the case of a fiat fi, spectrum galaxy with a complete break at the Lyman 
continuum edge, which should be representative of typical high redshift galaxies. For each of 
5 colors (Johnson B and V, Kron-Cousins R and / and Sloan z) we computed the convolution 
of the spectrum with the filter transmission and the forest transmission and computed the 
decrement relative to a flat fi, object without a break. For regions of the spectrum where the 
Lya forest would lie in the filter we interpolated the numbers in Table 2a to determine the 
average transmission, and similarly for the Ly/3 transmission from Table 2b. For wavelengths 
between the Lyman continuum break and the Ly7 hue we used the results of Table 2c. The 
results are relatively insensitive to this approximation of the higher order Lyman series 
scattering. The results are shown in Figure 4 and summarised in Table 3 where the changes 
in magnitudes produced by the decrement are given. The difference between the magnitude 
shifts in any two fllters gives the color of the galaxy at the redshift in the AB magnitude 
system of Oke. Therefore, a,t z — 5.7, for example, the AB {R — z) color is 2.19. Similar 
results may be computed for any input spectrum using the transmissions in the tables. 



5. Conclusions 

We have measured the transmission of the Lya forest produced by neutral hydrogen 
scattering in the intcrgalactic medium between redshifts 2 and 6.3 using high signal to noise, 
high resolution {R > 5000) observations of 50 quasars spread over the redshift range. The 
transmitted fractions show a relatively smooth evolution over the entire range of redshifts, 
which can be modelled with a smoothly decreasing ionization rate. We have used the tab- 
ulated transmissions to calculate the colors of a flat-spectrum galaxy with a large intrinsic 
Lyman continuum break. 
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Table la. High Redshift ESI Observations 

Quasar Mag. Expo, (hrs) Zem 

SDSS 0206+1216 19.9 3.0 4.81 

SDSS 0211-0009 20.0 7.0 4.90 

SDSS 0231-0728 19.2 4.33 5.42 

SDSS 0338+0021 20.0 7.25 5.01 

BR 0353-3820 18.0 2.5 4.58 

PSS 0747+4434 18.1 3.0 4.39? 

FIRST 0747+2739 17.2 1.5 4.17 

SDSS 0756+4104 20.1 3.0 5.09 

SDSS 0836+0054 18.8 6.0 5.82 

BRI 0952-0115 18.7 4.25 4.42 

SDSS 1030+0524 19.7 6.28 6.28 

BR 1033-0327 18.5 2.0 4.51 

SDSS 1044-0125 19.7 5.75 5.74 

SDSS 1048+4637 18.4 6.23 6.23 

SDSS 1148+5251 19.1 6.39 6.39 

SDSS 1204-0021 19.1 3.0 5.07 

BR 1202-0725 18.7 3.0 4.61 

SDSS 1306+0356 19.6 6.75 5.99 
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Table la — Continued 



Quasar Mag. Expo, (hrs) Ze.m 



SDSS 1321+0038 20.1 2.0 4.71 

SDSS 1605-0122 19.4 3.75 4.93 

WFSJ 1612+5255 19.9 1.5 4.95 

SDSS 1737+5828 19.3 7.5 4.85 

SDSS 2200+0017 19.1 5.3 4.78 

SDSS 2216+0013 20.3 3.0 5.00 

BR 2237-0607 18.3 11.6 4.55 
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Table lb. Low Redshift HIRES Observations 

Quasar Mag. Expo, (hrs) Zem 

Q0014+813 16.5 7.8 3.38 

HS0119+1432 16.7 7.8 2.87 

Q0256-000 18.7 5.5 3.37 

Q0302-003 18.4 10.7 3.28 

Q0636+680 16.5 7.4 3.18 

HS0741+4741 17.5 8.0 3.22 

HS0757+5218 17.5 4.0 3.2 

Q08279+5255 15.2 3.8 3.87 

Q0956+122 18.2 7.3 3.30 

Q1032+0414 18.1 4.0 3.39 

J1057+4555 17.7 11.3 4.10 

HE1104-1805 16.3 2.0 2.31 

HE1122-1648 16.6 3.5 2.40 

Q1159+123 17.5 8.0 3.50 

Q1206+119 17.9 5.0 3.113 

B1310+4254 18.4 5.2 2.561 

Q1422+2309 16.7 16.2 3.63 

Q1623+269 16.0 4.0 2.52 
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Table lb — Continued 



Quasar Mag. Expo, (hrs) Zem 



FBQS1634+3203 16.6 4.0 2.35 

HS1700+6416 16.1 11.0 2.748 

Q1937-1009 17.5 9.3 3.805 

HS1946+7658 15.8 12.0 3.02 

Q2000-330 17.6 4.7 3.77 

Q2126-158 17.3 4.4 3.28 

HE2347-4342 16.1 6.7 2.88 
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Table 2a. Mean Transmission Blueward of Lya 

^median (-^) 2^max -^min 

2.40 0.820 0.900 0.710 

2.46 0.809 0.879 0.735 

2.54 0.818 0.899 0.761 

2.61 0.778 0.833 0.700 

2.67 0.728 0.839 0.608 

2.72 0.800 0.893 0.702 

2.75 0.780 0.837 0.700 

2.79 0.764 0.811 0.648 

2.83 0.698 0.806 0.575 

2.86 0.695 0.806 0.552 

2.89 0.729 0.808 0.643 

2.92 0.694 0.799 0.492 

2.95 0.788 0.847 0.744 

2.98 0.712 0.812 0.644 

3.01 0.723 0.833 0.618 

3.03 0.692 0.821 0.456 

3.08 0.697 0.774 0.639 

3.09 0.719 0.855 0.545 
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Table 2a — Continued 



-2'median (-^) 2^max -Tmin 

3.14 0.767 0.849 0.711 

3.19 0.717 0.785 0.650 
3.24 0.675 0.713 0.630 
3.35 0.687 0.800 0.597 
3.48 0.628 0.730 0.532 
3.59 0.597 0.729 0.467 
3.65 0.553 0.633 0.449 
3.75 0.495 0.574 0.414 
3.84 0.495 0.624 0.352 
3.89 0.563 0.803 0.290 
3.94 0.487 0.612 0.363 
3.99 0.467 0.565 0.346 
4.03 0.388 0.469 0.302 
4.08 0.434 0.577 0.297 
4.12 0.491 0.640 0.367 
4.17 0.487 0.554 0.419 

4.20 0.441 0.577 0.342 
4.23 0.414 0.492 0.348 



Table 2a — Continued 



-2'median (-^) 2^max -Tmin 



4.26 0.404 0.487 0.309 

4.29 0.300 0.453 0.135 

4.31 0.439 0.588 0.326 

4.35 0.350 0.420 0.279 

4.37 0.398 0.551 0.219 

4.41 0.426 0.621 0.288 

4.44 0.348 0.397 0.253 

4.48 0.328 0.401 0.277 

4.52 0.316 0.408 0.187 

4.56 0.376 0.472 0.247 

4.59 0.297 0.366 0.148 

4.66 0.222 0.317 0.110 

4.68 0.215 0.273 0.093 

4.74 0.238 0.303 0.194 

4.82 0.246 0.334 0.200 

4.98 0.182 0.222 0.151 

5.13 0.147 0.321 0.095 

5.26 0.084 0.208 0.016 
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Table 2 a — Continued 



-^-median (-^) 2^max 2^min 

5.42 0.090 0.150 0.027 

5.51 0.060 0.161 0.025 

5.60 0.082 0.178 0.044 

5.70 0.028 0.071 -0.001 

5.87 0.009 0.025 -0.001 

6.05 0.0006 0.007 -0.002 
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Table 2b. Mean Transmission Blueward of Ly/3 



-2'median (-^) ^max 



mm 



3.06 
3.18 
3.27 
3.52 
4.05 
4.30 
4.42 
4.59 
4.71 
4.81 
4.95 
5.66 
6.11 



0.735 
0.718 
0.707 
0.641 
0.502 
0.502 
0.457 
0.378 
0.253 
0.243 
0.258 
0.082 
0.0002 



0.820 
0.785 
0.801 
0.700 
0.611 
0.582 
0.560 
0.625 
0.413 
0.328 
0.353 
0.140 
0.009 



0.627 
0.624 
0.616 
0.585 
0.394 
0.433 
0.380 
0.200 
-0.013 
0.085 
0.162 
0.035 
-0.008 
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Table 2c. Mean Transmission Blueward of Ly7 



-^-median (-^) 2^max 2^min 

3.32 0.715 0.749 0.650 

4.04 0.428 0.668 -0.422 

4.49 0.449 0.541 0.373 

4.78 0.249 0.549 -0.039 

5.00 0.240 0.416 0.173 

6.15 0.038 0.126 -0.001 
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Table 3. Colors of High redshift Galaxies and 
Quasars 



Redshift z I R V B 



3.0 0.00 0.00 

3.1 0.00 0.00 

3.2 0.00 0.00 

3.3 0.00 0.00 

3.4 0.00 0.00 

3.5 0.00 0.00 

3.6 0.00 0.00 

3.7 0.00 0.00 

3.8 0.00 0.00 

3.9 0.00 0.00 

4.0 0.00 0.00 

4.1 0.00 0.00 

4.2 0.00 0.00 

4.3 0.00 0.00 

4.4 0.00 0.00 

4.5 0.00 0.00 

4.6 0.00 0.00 

4.7 0.00 0.00 



0.00 0.00 0.37 

0.00 0.01 0.48 

0.00 0.04 0.59 

0.00 0.08 0.65 

0.00 0.13 0.68 

0.00 0.18 0.72 

0.00 0.25 0.84 

1.10 0.33 1.00 

0.00 0.43 1.17 

0.02 0.51 1.41 

0.07 0.56 1.76 

0.14 0.60 2.17 

0.21 0.65 2.68 

0.29 0.69 3.55 

0.39 0.72 5.39 

0.50 0.85 7.80 

0.62 1.01 9.19 

0.79 1.29 10.5 
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Table 3 — Continued 



Redshift z I R V B 



4.8 0.00 0.00 

4.9 0.00 0.00 

5.0 0.00 0.03 

5.1 0.00 0.10 

5.2 0.00 0.18 

5.3 0.00 0.29 

5.4 0.00 0.42 

5.5 0.00 0.56 

5.6 0.00 0.72 

5.7 0.00 0.91 

5.8 0.00 1.16 

5.9 0.01 1.46 

6.0 0.04 1.87 

6.1 0.12 2.57 

6.2 0.23 3.30 

6.3 0.35 3.66 

6.4 0.49 3.93 



0.99 1.63 12.2 

1.13 1.97 13.6 

1.19 2.30 18.2 

1.24 2.56 25.0 

1.30 2.92 25.0 

1.39 3.42 25.0 

1.51 4.06 25.0 
1.69 5.31 25.0 
1.94 6.98 25.0 
2.19 8.71 25.0 

2.52 10.0 25.0 
2.97 11.1 25.0 
3.38 12.3 25.0 
3.83 14.0 25.0 
4.34 20.4 25.0 
4.83 25.0 25.0 
5.45 25.0 25.0 
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Fig. la. — Transmitted flux blueward of Lya emission as a function of redshift for both the 
ESI and HIRES samples. In all cases, flux was computed in 8 bins of width 15 A between 
1080 A and 1185 A. For the z > 4.5 ESI sample (filled squares), the flux was normalized to 
a quasar power law continuum with slope —1.25 normalized at 1350 A. Large open squares 
denote data for real or suspected BAL QSOs in the sample, namely SDSS 1605—0112, 
SDSS 1044-0125, and SDSS 1048+4637. For the lower redshift HIRES sample (crosses), 
the flux was normalized to a continuum that was interpolated locally between transmission 
peaks. For the z < 4.5 ESI sample, the continuum was computed by both methods and 
is plotted (diamonds) assuming the local interpolation. See text for a discussion of the 
comparison. Where the transmitted flux is less than 0.0015, we have shown it at this nominal 
value. 
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Fig. lb. — As in Fig. la for transmission blueward of Ly/3. Flux was computed in two bins 
of width IsA between 980A and lOlOA. 
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Fig. 2a. — Optical depth blueward of Lya emission as a function of redshift for the com- 
plete sample of Figure la. Starting with the highest redshift points, mean transmissions 
were computed in bins of six data points. The optical depths corresponding to the mean 
transmissions are plotted (filled squares) as a function of the median redshift in each bin. 
The error bars show the range of extremes of transmission within each bin, converted to 
optical depth. Mean transmissions and median redshifts in each bin are tabulated in Table 
2a, along with the redshift range in each bin. The dashed line shows a uniform ionized 
Gunn-Peterson model, normalized to match the observations. The solid line is the model of 
Equation (2), with g set equal to 0.74 x ((1 + z)/6)~'^-^. Open symbols denote Lya optical 
depths from Figure 6 of White et al. (2003) for the quasars SDSS 1148+5251 (diamonds) 
and SDSS 1030+0524 (squares). 
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Fig. 2b. — As in Fig. aA for Ly /3. The dashed hne is a uniform ionized Gunn Peterson 
model for Ly /3, computed with the same normahzation as in Figure 2a. It underpredicts 
the Ly f3 absorption because of the inhomogeneous nature of the IGM. The sohd hne is the 
computed absorption from Equation 3 for the same run of ionization as in Figure 2a. See 
text for discussion. 
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Fig. 3. — Evolution of the normalized ionization rate as a function of redshift computed from 
mean transmissions using Equation 2 {filled squares). Mean transmissions and dispersions 
and median redshifts were computed in 12-point bins, starting at the highest redshift. The 
horizontal error bars show the range of redshift in each bin and the vertical error bars show 
the ionization parameter corresponding to ±1 a in the mean transmission in each bin. The 
open squares are the measurements of McDonald et al. (2000) at lower redshift. The dashed 
hne is a power law of the form 0.74 x ((1 + z)/6)~^-^, fitted to the present data for z > A. 
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Fig. 4. — Magnitude change expected as a result of IGM opacity as a function of redshift for 
a flat fiy spectrum galaxy with a complete break at the Lyman continuum edge, for each of 5 
colors: Johnson B and V, Kron-Cousins R and / and Sloan z. The Suprime-Cam filter trans- 
missions (http://www.noaj.org/Observing/Instruments/SCam/sensitivity.html) were used 
for the specific transmissions adopted. Approximate forest decrements were computed from 
the Lya, Ly/3, and Lyj transmissions of Table 2. 



